Different particle alignments in N ~ Z Ru isotopes studied by the shell model 
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Experimentally observed heaviest N « Z nuclei, Ru isotopes, are investigated by the shell model 
on a spherical basis with the extended P + QQ Hamiltonian. The energy levels of all the Ru 
isotopes can be explained by the shell model with a single set of force parameters. The calculations 
indicate enhancement of quadrupole correlations in the N — Z nucleus 88 Ru as compared with the 
other Ru isotopes, but the observed moments of inertia seem to require much more enhancement 
of quadrupole correlations in 88 Ru. It is discussed that the particle alignment takes place at 8 + in 
90 Ru but is delayed in 88 Ru till 16 + where the simultaneous alignments of proton and neutron pairs 
take place. The calculations present interesting predictions for 89 Ru that the ground state is the 
1/2" state and there are three A J = 2 bands with different particle alignments including the T — 
p — n pair alignment. 

PACS numbers: 21.10Hw,21.10.Re,21.60.Cs,23.20.Lv 



I. INTRODUCTION 

The so-called delay of alignment in the N = Z even- 
even nuclei is observed in the 64 < A < 88 region 0, 0, 
0, Q and in the lighter nucleus 48 Cr. This phenomenon 
is a sign of strong proton-neutron (p — n) correlations 
in the same shell |5|, and the special collectivity in the 
N = Z even-even nuclei suggests a strong collaboration 
of p — p, n — n and p — n correlations in the A = Am 
nuclei with N = Z = 2m which can be called the a- 
like (T = 0) 2p - 2n correlations 0,011. A theoretical 
investigation of this is challenging, which belongs to the 
study of the propert ies of the y — n interaction in N = Z 
nuclei [1 liH TTTl Til 111 111 Il5|. The experimental study 
of heavy N = Z nuclei has reached 88 Ru 3, 4]. The new 
data have revealed that there is a remarkable difference 
between neighboring even-even nuclei with N = Z and 
N = Z + 2 in the lgg/2-subshell region. The qualitative 
difference between 88 Ru and 90 Ru ( 84 Mo and 86 Mo) in 
the backbending plots of the yrast bands is different from 
the conditions in lighter nuclei Zr, Sr, etc., and casts a 
new light on the problem of the delayed alignment. 

A theoretical explanation of the delayed alignment for 
heavy N = Z nuclei 84 Mo, 88 Ru etc. is presented in Refs. 
0, El with the projected shell model on the deformed 
basis 0, 0| . The projected shell model reproduces 
the graphs of observed moments of inertia, by adopting 
commonly accepted deformations for those nuclei. The 
adopted deformations manifest that the deformation is 
larger for even-even N — Z nuclei when compared with 
N > Z nuclei. In other words, the delayed alignment in 
88 Ru is related to the large deformation. The study with 
the projected shell model 0, ^| suggested an enhance- 
ment of the p — n quadrupole-quadrupole (QQ) interac- 
tion in the N = Z nuclei. It is our interest to understand 
the structural difference between the N = Z even-even 



nucleus 88 Ru and neighboring isotopes in various aspects. 
In this paper, we make the study using the shell model 
calculations on the spherical basis which is free from fix- 
ing the deformation parameter. 

The extended P + QQ model 0, |2(j reproduces 
observed energy levels and B(E2) in N w Z I/7/2- 
subshell nuclei and is capable of describing the backbend- 
ing phenomena. It has successfully clarified characteris- 
tics of the structure of a heavier N — Z nucleus 64 Ge 
in the configuration space (2p 3/2 , 1/5/9, 2p 1/2 , 159/2) 8 
in a recent shell model calculation |2l|. The heav- 
iest N = Z nucleus experimentally observed, 88 Ru, 
which is expected to have the approximate configuration 
(2jt?3/2, 1/5/2, 2pi /2 , l59/2)~ 12 , is a good target to study 
the delayed alignment using the shell model calculation. 
The success in 64 Ge suggests that the extended P + QQ 
model provides a reliable interaction for the study of the 
heavy N sw Z nuclei. We carry out shell model calcula- 
tions using the extended P + QQ model with a single set 
of force parameters fixed for 88 Ru and heavier Ru iso- 
topes. The calculations, which are carried out with the 
calculation code |22j . have huge dimensions (maximum 
dimension is 165 x 10 6 for 88 Ru) and can be regarded as 
realistic ones. We investigate the structure of Ru isotopes 
and examine whether the difference between 88 Ru and 
90 Ru is reproduced or not by the spherical shell model, 
in Section III. The present shell model predicts interest- 
ing features of the odd- A isotope 89 Ru between 88 Ru and 
90 Ru. The prediction for 89 Ru is shown in Section IV. 

Since the delayed alignment in 88 Ru seems to be re- 
lated to the strong quadrupole correlations and the large 
quadrupole deformation [j, , we pay our attention to 
the roles of the QQ force which induces the quadrupole 
correlations and deformation. It is interesting to see the 
competition between the like-nucleon (p — p and n — n) 
interaction and p — n interaction of the QQ force. We 
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also examine a possible contribution of the isovector QQ 
force to the properties of Ru isotopes. 



II. THE MODEL HAMILTONIAN 



The extended P + QQ Hamiltonian is given by 
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where e a is a single-particle energy, H mc denotes the 
monopole corrections, Pjmtk. is the pair operator with 
angular momentum J and isospin T, and Q2M (O^m) 
is the isoscalar quadrupole (octupole) operator (see Ref. 
[20j). The force strengths X2 an d X3, are defined so as 
to have the dimension of energy. Following Ref. [2l| , we 
adopt the model space (2p 3 / 2 , I/5/2, ^Pi/2, 1^9/2) an d in- 
troduce the isoscalar octupole-octupole force Hoo'- Note 
that the Hamiltonian is isospin-invariant and includes the 
p — n pairing forces in addition to the p — n QQ force. 

The isoscalar QQ force Hqq can be divided into three 
parts, p — p, n — n and p — n parts. We shall use the 
notations x% P , X%nn an d X%n f° r their force strengths. 
In terms of the pair operators Pjmtk, the isoscalar QQ 
force is expressed as 

9 E ~kt E W JT=\ P JMlK P JMl K 

,1 moo, (2) 
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where xz = 



X°>- The symbol Wjt proportional 
to the Racah coefficient really has four subscripts related 
to the four orbits of c^c^csc^. The first line of Eq. © 
brings about the isovector pairing interactions, where the 
strengths x\ =1 , x\__ 1 and x}. =Q stand for the n — n, 
p — p and p — n interactions. The second line of Eq. 
(|2~|) brings about the isoscalar p — n pairing interactions. 
The p — n part of Hqq is enhanced by enlarging the 

force strength x% n ( x \=o an d 2^=0) i n Refs. 0,0 EI- 
The Hamiltonian ceases to be isospin-invariant, with the 
isospin not being a good quantum number there. 

There is a possibility of the isovector QQ force con- 
tributing to the collective motion in the heavy N « Z 
nuclei. The isovector QQ force Hqq" 1 with the force 
strength x\ is & lso rewritten in the same form as Eq. 

-3x1 We 



can write the sum of Hqq and Hqq 1 in the same form 



(0) with the relations x\ — x\ an d a;°=o 

ite the sum of 
as Eq. 10, where x\ 
3nsk 

the following combination of the interaction strengths 

(l + a)x, xl = -ax, (3) 



X2 + X2 and x ° K =o = X2- 3 X2- 
If we consider a restricted sum of Hqq and Hq^ 1 with 
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the QQ force is written as 
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By changing the mixing parameter a, we can enhance 
the p — n part of the QQ force, which corresponds to the 
isoscalar p — n pairing interactions in the second line of 
Eq. without violating the isospin invariance of the 
Hamiltonian. 



III. DIFFERENCE BETWEEN 88 RU AND 90 RU 

Using the extended P + QQ Hamiltonian we 
carried out shell model calculations in the hole s pac e 
(i#s>/2> 2Pi/ 2 ' I/5/2' ^3/2) w ^ n tne calculation code |22| . 
The single-hole energies e\ depend on H mc and the force 
strengths as well as e a through the hole transformation. 
We treated the hole energies ejj as parameters instead of 
the single-particle energies e Q . We tried various combina- 
tions of the parameters e^, H mc , go, 92, X2 an d X3; an d 
determined these parameters so as to reproduce overall 
energy levels of the Ru isotopes. The adopted parameters 
are 



fc 9/2 



0.0, 



Jl 

■'1/2 



1.1, 



£ 5/2 



5.5, 



£ 3/2 



6.0, 



.go = 0.26(92/A), g 2 = 0.12(92/A) 5 / 3 , 



18- 
17 



16- 
16 

15" 

14 ; 

5\r 13 



4r 



n 3 



92i-> 

Ru 

-18 18" 

-17 17- 

"16 

"15 

16" 



94 



Ru 



12- 



10- 



-13 

-11 
-9 
"7 
-5 



15 - 
14 - 
13 



12- 
10 - 



-17 

-16 
-15 



-13 

-11 
-9 

-7 
-5 



15 ■ 

14 ■ 
13 ■ 



12 " 
10 ■ 



-15 



-13 



11 12 
9 

v 10- 



-13 



-1 1 

"9 



exp 



exp 



FIG. 1: Energy levels of 92 Ru and 94 Ru. The label "A" stands 
for the energy levels calculated with the parameter set (5) and 
"exp" for the observed ones. 
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xl 



0.26(92/A) 5 / 3 , xl = 0.04(92/A) 2 in MeV, (5) 



and i? mc is fixed at zero (the J-independent isoscalar 
monopole term is not determined, because we do not deal 
with the binding energy in this paper). Changing the 
monopole corrections H mc does not significantly improve 
the energy levels. The relative position of e^ 2 and e\, 2 
is responsible to that of the positive and negative parity 
states. In our trials, the values of £g , 2 and £3 , 2 listed 
in (J5J are best and the exchange of the two values does 
not improve the energy levels. The hole levels l/g/ 2 and 

2py 2 seem to lie far from 2p'y 2 . This is the reason why 
the subspace (2p 1 / 2 , I59/2) works well for A > 86 nuclei 
in Refs. |2^,[24|. The force strengths g , g 2 , X2 an d X3 
adopted are similar to those used in the study of 64 Ge 

m 

The parameter set (0 reproduces well the energy lev- 
els (the patterns and order of the positive- and negative- 
parity levels) of Ru isotopes, not only the even- A nuclei 
92 Ru and 94 Ru but also the odd- A nuclei 91 Ru and 93 Ru 
as shown in Figs. ^andE! The agreement between the- 
ory and experiment for the odd-parity state is worse than 
that for the even-parity states. The calculation, how- 
ever, reproduces the observed energies within the error 
0.8 MeV. 
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FIG. 2: Calculated and observed energy levels of 91 Ru and 
93 Ru. The spin of each state is denoted by the double number 
2J. 



A. Dependence on the QQ force strengths 
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FIG. 3: Comparison of calculated energy levels with observed 
ones for 88 Ru. The calculated results are obtained with the 
different strengths of the QQ force, A, B, C, D and E. 



the difference between 88 Ru and 90 Ru in the backbending 
plot (we call it "J — uj graph") as shown in Figs. Eland 
[S] The calculation reproduces the sharp backbending 
at J = 2j — 1 = 8 observed in 90 Ru and shows no clear 
backbending in low-spin states of 88 Ru in agreement with 
the experiment. 

The backbending plots, however, reveal insufficiency 
for the most collective low-lying states. The results A do 
not reproduce the slope of the J—lu graph up to J = 8 for 
88 Ru and up to J = 6 for 90 Ru. The slopes of the J — u> 
graphs for the collective bands are considerably affected 
by the strength X2 °f t ne QQ force Hqq above all other 
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The energy levels obtained for 88 Ru and 90 Ru, which 
are shown in the column A of Figs. |3|and21 are consistent 
with the observed ones. The parameter set A describes 



FIG. 4: Comparison of calculated energy levels with observed 
ones for 90 Ru. The calculated results are obtained with the 
different strengths of the QQ force, A, B, C, D and E. 
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FIG. 5: The J — u graph of Fig. 1. 



FIG. 6: The J — u graph of Fig. 2. 



force strengths. This is naturally understood, because 
the moment of inertia of a rotational band depends on 
the magnitude of deformation and the QQ force drives 
the quadrupole deformation. Let us try to improve the 
J — u> graphs for 88 Ru and 90 Ru by readjusting the QQ 
force strength. 

We first strengthen the p — n QQ interaction by adding 
the isovector QQ force Hq^ 1 in the form Q so as to con- 
serve the isospin invariance. The results obtained with 
the mixing parameter a = 0.125 (see Eq. ©) are shown 
by the notation B in Figs. - El The J — ui graph is 
improved for 90 Ru. For 88 Ru, the result B removes the 
slight backbending at J = 8 of the result A. The param- 
eter set B reproduces quite well the overall energy levels 
of the Ru isotopes 90 Ru, 91 Ru, 92 Ru, 93 Ru and 94 Ru. 
For the high-spin states, however, the parameter set A 
is better than B. The change from A to B pushes up 
the high-spin levels higher as the spin J increases. Since 
the configuration (l3g/ 2 ) m is dominant in the high-spin 
states, the inadequacy for the high-spin states suggests 
that the enhanced p — n QQ force strength {x^ =0 = 1.5x) 
is too strong for the lgg/2 subshell. The remaining devi- 
ation of the calculated J — u) graph from the experimen- 
tal one for 88 Ru indicates room for improvement in the 
model space and in the interactions of our model. 

Results similar to those of B are obtained by strength- 
ening the p — n part of the isoscalar QQ force Hqq 

(by enlarging x% n )- Th e results of xlvn = l- 25 X° P p are 
shown by the notation C in Figs. |3]-|6] (Note that the 
p — n force strength x%n = ^-^X2p P is used by Sun et 
al. to increase the particle alignment frequency for 
88 Ru and 84 Mo.) Although the high-spin levels of C are 
pushed up a little higher as compared with those of B, 
the parameter sets B and C yield similar results, not only 
for the energy levels but also for B(E2) values and the 
quadrupole moment Q with respect to the yrast bands 
of 88 Ru and 90 Ru as shown later on. We also tried to 



strengthen all of the p — p, n — n and p — n QQ interac- 
tions. The enlargement of X2 to 1.1 x X2 yields results 
similar to those of B and C. The results are denoted by D 
in Figs. [31- EI Within the small increase of the QQ force, 
we have found no evidence that the p — n QQ interaction 
is stronger than the p — p and n — n QQ interactions, 
and there is no choice between the isospin-variant and 
-invariant enhancements of the p — n QQ interaction, as 
well. 

It should be noted here that the strength of the p — n 
interaction does not directly correspond to the strength 
of the p — n correlations. According to the single- j shell 
calculation with the extended P+QQ force [25j , the p—n 
correlation energy becomes largest at N = Z in nuclei 
with the same Z even though the same p—n interaction is 
used for those nuclei, while the p— p and n — n correlation 
energy does not show such a specific feature. 

In Fig. the discrepancy between the calculated mo- 
ments of inertia (B, C and D) and observed ones is still 
large for 88 Ru. The calculations B, C and D cannot suf- 
ficiently reproduce the observed large angular frequency 
at J = 8. If we want to obtain a better slope of the J — lu 
graph for 88 Ru, we must enhance the QQ force strength 
X2 much more. The slope of the J — u> graph observed 
in 88 Ru cannot be well reproduced even by strengthen- 
ing the p — n QQ interaction further in the way B or 
C. Results obtained with the isoscalar QQ force Hqq 
strengthened by 1.5 (1.5 x xV) are denoted by E in Figs. 
13 - El The calculation E reproduces the large angular 
frequency and the slope (moment of inertia) of the J — uj 
graph for 88 Ru but yields rather bad results for 90 Ru. 
This suggests that the collectivity of the quadrupole cor- 
relations is different between 88 Ru and 90 Ru. The ob- 
served moment of inertia and large angular frequency 
indicate a stable rotation of 88 Ru, while the energy lev- 
els and the sharp backbending at 8 + reveal a deviation 
from the rotation in 90 Ru. Our spherical shell model 
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calculation predicts that clear backbending does not oc- 
cur up to J = 14 in 88 Ru, as predicted by the projected 
shell model 0, 0] . The present results indicate a spe- 
cial enhancement of the quadrupole correlations in 88 Ru 
in contrast to the other Ru isotopes with N > Z. This 
is consistent with the results obtained by the projected 
shell model 0, Il6| , in which the deformation parameter 
is fixed to be 0.23 for 88 Ru and 0.16 for 90 Ru (the former 
is 1.4 times as large as the latter). Our spherical shell 
model requires the enhancement of the QQ force instead 
of the enlargement of the deformation for the N = Z nu- 
cleus 88 Ru. This suggests that the present model space 
(2p 3 / 2 , 1/5/2, 2pi/ 2 , 139/2) may not be sufficient and the 
lower orbit I/7/2 or the upper one 2d 5 / 2 should be in- 
cluded possibly. 



B. Difference between 88 Ru and 90 Ru in structure 

The backbending at J = 2j — 1 = 8 in the I39/2- 
subshell nucleus 90 Ru is contrast to no backbending at 
J — 2j — 1 — 6 in the l/ 7 / 2 -subshell nucleus 50 Cr, while 
the resistance to the backbending is common to the N — 
Z nuclei 88 Ru and 48 Cr. Let us discuss the difference 
between 88 Ru and 90 Ru in the structure which appears 
in the J — ui graphs of Figs. [S]and[fjl In Tables HI and ITU 
we tabulate the expectation values of proton and neutron 
numbers (n a ) in the respective orbits for the yrast states 



TABLE I: Expectation values of proton (neutron) numbers in 
the respective orbits for the yrast states of 88 Ru calculated 



of 88 Ru and 90 Ru. The tables show that more protons 
jump up from the pf subshell to the l<?g/ 2 one in 88 Ru 
than in 90 Ru, and the same is true for neutrons if the 
extra neutron pair is subtracted from the neutron number 
(n g9 / 2 ) for 90 Ru. Two things are characteristic in the 
N w Z Ru isotopes, which is different from the situation 
of the N w Z Cr isotopes in the I/7/2 subshell. First, 
the l<79/2 subshell where the Fermi level lies is just above 
the pf subshell and there is a considerably large degree 
of freedom for 139/2- Secondly, the two subshells have 
opposite parities. These conditions permit only nucleon 
pairs jumping up to I39/2 and induce strong p — p, n — n 
and p — n correlations in lg 9 / 2 . We can suppose that the 
collaboration of the p — p, n — n and p — n correlations 
results in the a-like 2p — 2n correlations especially in the 
139/2 subshell, in the N = Z nucleus 88 Ru where the 
p — n correlations are enhanced. 

The sharp backbending at 8 + in Fig. for 90 Ru coin- 
cides with the increase of (n g g/2) and decrease of (71.pi/2) 
for neutron (and their decrease and increase for proton) 
at J = 8 in Table [H] This change is explained by the 
alignment of a neutron pair in 139/2- The results A and 
B in Table [H] suggest the following explanation: There 
is an extra neutron pair which cannot form a T = 2p- 
2n quartet |3, |g, |26j in 90 Ru. The extra neutron pair 
has a dominant probability to be a pair with J — and 
T = 1 , and contributes to the collective 2p — 2n corre- 



TABLE II: Expectation values of proton and neutron numbers 
in the respective orbits for the yrast states of 90 Ru calculated 
with the different strengths of the QQ force, A and B. 

cal. A 



strengths of the QQ force, A and B. 
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20 
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lations through the exchange with a neutron pair in the 
quartets. The excitation till 6 + owes to the motion of the 
quartets. At 8 + , the extra neutron pair aligns the angular 
momentum to be J = 9/2 + 7/2 ( J = 8, T = 1) in lg 9/2 
and breaks away from the collective 2p — 2n correlations, 
which increases the 1 179/2 neutron number. The weak- 
ened 2p — 2n correlations somewhat hinder proton pairs 
jumping up to the lgg/2 subshell from the pf sub-shell, 
which decreases the lc/9/2 proton number. The result A 
for 88 Ru in Table [I] shows a similar sign at J = 8, but 
the observed J — u> graph denies such a pair alignment 
in 88 Ru. The calculated results B, C, D and E which 
are better for the very collective low-lying states sweep 
away the sign of a structural change in (n a ) (the result 
B is shown in the lower part of Table . By combining 
our result for 88 Ru with the J — uj graphs experimentally 
observed in other N = Z nuclei, we can say that the one- 
pair alignment is hindered in the N = Z nuclei due to the 
strong 2p — 2n correlations. This may be the reason for 
the durable increase of angular frequency in the N = Z 
nuclei. 

Instead, Fig. [31 in which the monotonous slope after 
J = 16 stands out, suggests a structure change at 16 + 
in 88 Ru. The projected shell model 0, E3 also predicts 
a backbending at J = 16 for 88 Ru. Moreover, the cal- 
culated result A in Table [I] shows the increase of (n/5/2) 
and (ripi/2), and the decrease of (rigg/2) at J = 16. The 
same sign remains slightly in (n g9 / 2 ) of the result B and 
the sign disappears for the strong QQ force E. As men- 
tioned above, however, the enhanced QQ force of B, C, D 
and E is more or less too strong for the high-spin states. 
We can expect that the structural change at J = 16 will 
be observed in 88 Ru. This structural change seems to 
be caused by the simultaneous alignments of proton and 
neutron pairs at J = 2 x (9/2 + 7/2), since the strong 
2p — 2n correlations resist the single align ment of proton 
or neutron pair. The analysis in Ref. |27| . which predicts 
the simultaneous alignments of proton and neutron pairs 
at J = 2 x (7/2 + 5/2) without backbending due to the 
one-pair alignment in the I/7/2 N = Z nucleus 48 Cr, sup- 
ports our conjecture for the lgg/2 N = Z nucleus 88 Ru. 
This conjecture is also supported by the backbending to- 
ward J = 16 observed in 90 Ru. In Table ITT1 the increase 
of proton number (ng g / 2 ) &t J — 16 in the results A and 
B for 90 Ru suggests the proton-pair alignment in 1*79/2 
in addition to the neutron-pair alignment at J = 8. 



Effect of the 2d 5 / 2 orbit 



The large deformation on the deformed basis can be 
interpreted by the mixing of a large number of spherical 
single-particle orbits. The expansion of the configura- 
tion space instead of the enhancement of the QQ force 
is effective in our spherical shell model. The 2^5/2 orbit 
could contribute to the quadrupole correlations, because 
it strongly couples with the I39/2 orbit which plays a 
leading role in the Ru isotopes, through the large ma- 



trix element (I59/2IIQI 12^5/2). Adding the 2d 5 / 2 orbit 
to the model space (2p 3 / 2 , 1/5/2, 2pi/ 2 , I59/2), unfortu- 
nately, makes the number of the shell model basis states 
too huge. Instead of this, let us examine the contri- 
bution of the 2d 5 / 2 orbit within the truncated space 
(2pi/2j lfl'9/2) 2^5/2)- This space does not cause the spu- 
rious motion of the center of mass, and is expected to 
work well as the truncated space (1/7/2,2^3/2) without 
(2pi/2, 1/5/2) can explain the main features of the fr/2- 
subshell nuclei [28l [2(| because of the large matrix el- 
ement (I/7/2IIQH2P3/2) (the Q matrix element becomes 
large when Al = Aj = 2). We carried out the shell model 
calculations using the single-particle energies £1/2 = 0.0, 
£ 9/2 = 1-0 and £5/2 = 6.0 in MeV. The inclusion of 2d 5 / 2 
allows us to use weaker force strengths than those in 
Eq. ||SJ|. We replaced the A-dependence (92/A) x with 
(88/A) x for 50, 52, X2 and X3 in E q- ©■ The results for 
88 Ru and 90 Ru are shown by the dash lines (X) in Fig. 
□ 

In Fig. [7| the calculated J — w graph agrees well with 
the experimental one observed for the low-lying collec- 
tive states of 88 Ru. The agreement is better than those 
of B, C and D in Fig. (It is notable that Fig. [7]also 
predicts the alignment at J — 16 in 88 Ru.) This suggests 
that adding the 2d 5/2 orbit to (2p 3/2 , l/ 5 /2,2pi/ 2 , 139/2) 
with adjusted force strengths probably improves the cal- 
culated results for 88 Ru. The collectivity in the expanded 
space could enhance the quadrupole correlations without 
strengthening the QQ force. Since the large deformation 
causes a large admixture of the 2d 5 / 2 spherical orbit in 
the deformed Nilsson basis states, the present result is 
consistent with the prediction of Ref. 0, 0] that the 
deformation of the N = Z nucleus 88 Ru is large (0.23). 
In contrast to this, the inclusion of 2d 5 / 2 ruins the J — uj 



(Pl/2>99/2^V2) 



J (Pl/2>99/2> d 5/2) 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.1 

(Ej-Ej_ 2 )/2 (MeV) (Ej-Ej_ 2 )/2 (MeV) 



FIG. 7: The J — uj graphs obtained in the model space 
(2pi/2, l<?9/2, 2^5/2) for 88 Ru and 90 Ru, compared with the 
experimental ones. 
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graph for 90 Ru as shown in Fig. The energy levels 
obtained do not display the sharp backbending at J = 8 
but look like a stable rotation. The discrepancy says that 
the 2d 5 / 2 orbit must not so much join in the quadrupole 
correlations (or the 2d 5 / 2 orbit must be far from 1^9/2) 
and hence the deformation is not large for 90 Ru, which 
is consistent with a smaller deformation (0.16) adopted 
inRef. HG3. 

As mentioned in Ref. the QQ force gives in- 

verse magnitudes to the interaction matrix elements 

((59/2)i T =il^QQ l(59/2)j,T=i) with J = 6 and J = 8 
contrary to those of the ordinary effective interaction 
pflj . This defect has a bad influence on the J — u> graph 
at J = 8. If we replace the J = 6 and J = 8 matrix 
elements with those of Ref. [23| , the slight backbending 
at J = 8 in the calculated result for 88 Ru disappears 
as shown by the dot line (Y) in Fig. [7| Figures to 
are not free from the same influence, either. This defect, 
however, does not change the general situation. 



D. B[E2) and Q moment 

We have discussed the energy levels, J — u> graph 
and (n a ) so far. The B(E2) value and Q moment are 
good physical quantities to see the characteristics of the 
quadrupole correlations. We calculated the B(E2) values 
and Q moments for the yrast states of 88 Ru and 90 Ru, 
using the different strengths of the QQ force A, B, C, D 



TABLE III: B(E2 : Ji -> J/) for the yrast states of 88 Ru 
and 90 Ru calculated with the different strengths of the QQ 
force, A, B, C, D and E. The last column X shows the values 
obtained in the model space (2p x / 2 , I59/2, 2d 6 / 2 ). 





B(S2 : Ji -» 


Jf) 


in e z 


fm 4 


Ji — ► Jf 
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88 Ru 














2^0 


460 
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536 


522 


612 


510 


4^2 


630 


730 


722 


704 


832 


730 


6^4 
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780 


758 
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582 


8^6 
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833 


806 


964 


601 


10 8 


771 


897 884 862 


999 


619 


12 -> 10 


775 


896 


882 


859 


999 


568 


14 -> 12 


748 


880 


866 


844 982 


532 


16 -> 14 


671 


853 


839 


814 


962 


489 


18 -> 16 


622 


820 


806 


782 


931 


430 


20 -> 18 


71 


760 


748 


725 


872 


348 


yu Ru 














2^0 


296 


339 


348 


339 


482 


567 


4^2 


385 


451 


466 


456 


663 


793 


6^4 


294 


388 


418 


418 


698 


713 


8^6 


235 


236 


213 


234 


654 


657 


10 -> 8 


309 


322 


325 


331 


764 


673 


12 -» 10 


284 


290 


289 


270 


501 


621 


14 -> 12 


257 265 


227 218 


600 


538 


16 -> 14 


238 


274 


296 


299 


396 


440 


18 -> 16 


252 


275 


296 


300 


597 


413 


20 -> 18 


191 


199 


202 


205 


542 


335 



and E (corresponding to those in Figs. 13-0) and the 
model X (corresponding to X in Fig. 01. We used the 
effective charges e p = 1.5e and e„ = 0.5e, to compare 
the relative values of electric quadrupole quantities ob- 
tained with the different strengths of the QQ force. The 
calculated results arc tabulated in Tables ITTTI and llVl 

In Table UTTI the calculated B(E2) values of 88 Ru are 
much larger than those of 90 Ru, providing that the en- 
ergy levels of both nuclei are approximately reproduced. 
The ratios of the B(E2) values are more than 1.6 in the 
calculation A. In other words, the quadrupole correla- 
tions are much more enhanced in the N = Z nucleus 
88 Ru than in 90 Ru. This is consistent with the fact that 
a larger deformation is employed for 88 Ru as compared 
with 90 Ru in the projected shell model |3.IT^|. 

The modifications of the QQ interaction, B, C and D, 
somewhat enlarge the B(E2) values both for 88 Ru and 
90 Ru. The ratios of the B{E2) values for 88 Ru to those 
for 90 Ru are still large. The very strengthened QQ force 
E, which is required to reproduce the slope of the J — u 
graph for 88 Ru, enlarges the B(E2) values fairly for 88 Ru 
and drastically for 90 Ru. We have already seen that the 
strengthened QQ force E ruins the pattern of energy lev- 
els for 90 Ru. The enhanced B(E2) values in the column 
E of Table lTTTI are therefore too large for 90 Ru. We do not 
adopt the large B(E2) values in the column X for 90 Ru 
from the same reason, either. The quadrupole correla- 
tions must not be enhanced too much and the contribu- 
tion of the 2d 5 / 2 orbit should be small for 90 Ru. Namely, 
90 Ru may not be largely deformed. The truncated con- 
figuration space (2pi/ 2 , 159/272^5/2) yields B(E2) values 
comparable to those of the result A for 88 Ru, in spite of 
the small space. The expansion of the model space by 



TABLE IV: Quadrupole moment Q(J) 
states of 88 Ru and 90 Ru calculated with 
of the QQ force, A, B, C, D and E. 



in e-fm for the yrast 
the different strengths 



J 


A 


B 


C 


D 


E 


88 Ru 2 
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43.6 
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-22.5 
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16 


1.8 


4.1 


7.4 


8.0 


36.2 


18 


-0.7 
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33.5 


20 
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-3.9 


-2.7 


-2.0 


35.7 
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adding 2d 5/2 to (2p 3 / 2 , 1/5/2, 2pi/ 2 , 159/2) can make the 
B(E2) values larger, which could be appropriate to the 
enhanced quadrupole correlations in 88 Ru. 

The calculated quadrupole moments Q(J) tabulated 
111 Table show the same results as the B{E2) val- 
ues. From Table I1VI we can say as follows: The small 
enhancements of the p — n QQ interaction (B and C) 
change insignificantly the structure of 88 Ru and 90 Ru, 
while the strong enhancement of the QQ force by 1.5 
times (E) changes the structure of 88 Ru drastically. The 
large and roughly constant Q moments of 88 Ru suggest 
the quadrupole deformation. If the result E should not 
be adopted for 90 Ru, Table llVl and the energy levels insist 
that 90 Ru does not have a large deformation. 

The calculated B(E2) values and Q moments in Tables 
IIII and llVl testify the structural change due to the par- 



ticle pair alignment at 8 + in 90 Ru, in contrast to 88 Ru. 
The B(E2 : 8 + — > 6 + ) value decreases and Q moment in- 
creases at 8 + in 90 Ru, while the two values do not show 
any abrupt changes at 8 + in 88 Ru. On the other hand, 
the simultaneous alignments of proton and neutron pairs 
at 16 + leave a sign in the B(E2) values and Q moment 
in the result A both for 88 Ru and 90 Ru. 



IV. PREDICTION FOR 89 RU 

The 89 Ru isotope between 88 Ru and 90 Ru has not ex- 
perimentally been observed yet. Our model, however, 
predicts interesting features of 89 Ru. Figure |H1 shows the 
energy levels and relative B(E2) values obtained using 
the parameter set A for 89 Ru. The collective states con- 
nected by large B{E2) values are divided into four bands. 
They are the yrast states except for 15/2~ and 31/2+. 




11/2+ 



Exceptionally, we select the second states (15/2),/ and 
(31/2) + as collective states based on the B{E2) values 
and Q moments, which are adjacent to the yrast states 
(15/2)7/ and (31/2)+ respectively. The relative B(E2) 
values are denoted by the widths of the arrows in Fig. 
|H1 The inter-band E2 transitions which are not shown in 
Fig. |H1 are weak. 

It is remarkable that the predicted ground state of 
the middle l<7g/ 2 -subshell nucleus 89 Ru is the 1/2" state. 
This extraordinary event is reasonable from the system- 
atic lowering of the l/2~ state with decreasing N in 
odd- A Ru isotopes as seen in Fig. [21 Our model re- 
produces the systematic behavior of 1/2 - . Look at the 



TABLE V: Expectation values 
in the respective orbits for the 
with the QQ force strength A. 
tabulated. 



of proton and neutron numbers 
yrast states of 89 Ru calculated 
Calculated Q moments are also 



FIG. 8: Energy levels predicted for Ru. The widths of the 
arrows show the relative B(E2) values. 
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expectation values of proton and neutron numbers (n a ) 
for the band-head states 1/2", 3/2", 9/2+ and 7/2+ of 
89 Ru which are tabulated in Table This table shows 
that the 1/2" state has more protons in lgg/2 than the 
9/2+ state. From the comparison of Table|V|with Tables 
HI and ITT1 the l/2~ state resembles the ground state of 
88 Ru and the 9/2 + state resembles the ground state of 
90 Ru. Roughly speaking, the l/2~ state is constructed 
by adding one neutron to 88 Ru(0 + ), and the 9/2 + state 
by removing one neutron from 90 Ru(0 + ). The B(E2) val- 
ues of the negative parity bands larger than those of the 
positive parity bands indicate the stronger collectivity of 
the negative parity bands. This corresponds to the re- 
sult shown in Table|mlthat the B(E2) values of 88 Ru are 
larger than those of 90 Ru. From these comparisons, the 
difference between the 1/2" state and the 9/2 + state can 
be understood in terms of the a-like 2p — 2n correlations 
mentioned in the interpretation of the difference between 
88 Ru and 90 Ru. The strong a-like 2p — 2n correlations 
pull up more protons to lgg/2 in 1/2" than in 9/2 + , be- 
cause the disturbing extra neutron is absent in lgg/2 for 
the 1/2" state. The inversion of 9/2 + and 1/2" says 
that the a-like 2p — 2n correlations give a larger energy 
gain to the l/2~ state and the larger correlation energy 
compensates the energy loss of more nucleon jumps to 
the I39/2 subshell in 1/2" . 

We can expect that the A J = 2 bands on the 1/2 - and 
9/2+ states are similar to the ground-state bands of 88 Ru 
and 90 Ru, respectively. In fact, the J — lj graphs for the 
two bands Nl and PI of 89 Ru which are shown in Figs. 
151 and ITU1 are similar to those of 88 Ru and 90 Ru in Figs. 
and Figure 01 suggests no backbending at 17/2" 
(1/2" + 8) in the negative parity band Nl, while Fig. ^3 
predicts a backbending phenomenon at 25/2 + (9/2 + + 
8) in the positive parity band PI. The backbending at 
25/2+ in the band PI seems to be caused by the proton 
pair alignment parallel to the spin of the last odd neutron 
in l<7g/2, corresponding to the neutron pair alignment at 
8 + in 90 Ru. The increase of the proton number (n g9 / 2 ) at 
25/2+ testifies the proton pair alignment in 139/2- The 
small B(E2) value from 25/2+ to 21/2+ and the decrease 
of the calculated Q moment at 25/2+ show the structural 
change there. The coincident increase of the neutron 
number {n g g/ 2 ) at 25/2 + gives another evidence of strong 
p — n correlations in l<?g/2- For the negative parity band 
Nl, the value of B(E2 : 17/2- -> 13/2 - ) does not show 
any sign of such a structural change and the expectation 
values of proton and neutron numbers show no abrupt 
change, which corresponds to no backbending at 8 + in 
88 Ru. 

Figure however, predicts backbending at 33/2" 
(1/2" + 16) in the band Nl. The simultaneous increases 
of proton and neutron numbers (n s a/2) at 33/2" (see 
Table |Vj say that this backbending is due to the simul- 
taneous alignments of proton and neutron pairs (J = 16) 
in I39/2, corresponding to the four nucleon alignment in 
88 Ru. The small value of B(E2 : 33/2" -> 29/2") and 
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FIG. 9: The J — u> graph for the negative parity bands Nl 
and N2 of 89 Ru. The labels A and C stand for the parameter 
sets A and C. 
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FIG. 10: The J — u) graph for the positive parity band PI of 
89 Ru. 



the increase of the Q moment at 33/2" testify the struc- 
ture change. Figure shows a sign of another back- 
bending at 37/2 + in the band PI, that is possibly the 
alignment of two protons and three neutrons in 1 179/2, 
corresponding to the alignment at 16 + in 90 Ru. The spin 
of the last odd neutron is parallel to the spin of rotation 
or alignment in the band PI. 

Figure |5] shows also the J — uj graph for the negative 
parity band N2. This figure with Fig. [Hlsays that the low- 
lying collective states 3/2", 7/2", 11/2" and (15/2)^ of 
the band N2 and 5/2", 9/2", 13/2" and 17/2" of the 
band Nl are the partners in the angular momentum cou- 
pling 1/2" ® J (= J± 1/2). The similar B{E2) values 
and similar Q moments support the picture. The most 
remarkable backbending in 89 Ru takes place at the 19/2" 
state of the negative parity band N2. The small B(E2) 
values from 19/2" to (15/2)7/, (15/2)^ and 17/2" indi- 
cate a clear structure change at 19/2 - . The Q moment 
decreases abruptly at 19/2 - . For the J > 19/2 - states 
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of the band N2, the Q moments are nearly constant and 
so is the slope of J — ui graph. This phenomenon can- 
not be explained by the nucleon pair alignment coupled 
to J = 8 (T = 1), because the J = 19/2 state cannot 
be constructed by the coupling 1/2 <S> 8. If the structure 
change is due to a kind of alignment, the phenomenon is 
attributed to the p — n alignment J = 9 (T = 0). 

In 89 Ru, the efficient way to construct the J = 9 p — n 
pair is one proton jump to lg 9 /2- When one p — n pair 
aligns to J = 9 (T = 0) in I59/2, another pair which 
breaks away from the a-like 2p — 2n correlations is still 
possible to join in the monopole (J = 0, T = 1) pair- 
ing correlations, and to couple with the last odd nucleon 
in 2pi/ 2 to the total isospin T = 1/2. The decreases 
of neutron and proton numbers (^9/2) at 19/2 - testify 
the decline of the a-like 2p — 2n correlations due to the 
breaking away of the J = 9 p — n pair from a T = 
2p — 2n quartet. It should be noted that the disunion 
of a T = 2p - 2n quartet to the T — Q and T = 1 
pairs is prohibited for even-even nuclei. The T = p — n 
alignment at 19/2 - seems to be a unique phenomenon in 
the l(?g /2-subshell odd- A nuclei with N = Z ± 1 such as 
89 Ru. (The T = p — n alignment could take place in 
N = Z odd-odd nuclei.) In this connection, the small E2 
values from 21/2~ to 17/2~ is notable. It suggests that 
the p — n alignment ( J = 9, T = 0) contributes to the 
J > 21/2 - states of the band Nl. Actually, the states 
21/2- , 25/2" and 29/2" of the band Nl resemble the 
states 23/2~, 27/2~ and 31/2" of the band N2, with re- 
spect to the energy levels, expectation values of nucleon 
numbers (n a ), B(E2) values and Q moments. They could 
be members of the collective excitations coupled with the 
three nucleons 2pl /2 (lg% /2 lg% /2 ) J=9tT=a . 

The positive parity band P2 shows a rather compli- 
cated behavior. The very low- lying 7/2+ state is appar- 
ently related to the state of three nucleons with J = j — 1 
in a high-spin orbit j The small B(E2 : 19/2+ -> 

15/2+) value and the abrupt decrease of the Q moment 
at 19/2+ testify a structure change at the 19/2+ state. 

V. CONCLUSIONS 

We have carried out the shell model calculations on 
the spherical basis using the extended P + QQ Hamilto- 
nian with a single set of parameters in the model space 
(2^3/2, 1/5/2J 2pi/2, 139/2)- The calculations reproduce 
qualitatively well the overall energy levels observed in 
the Ru isotopes, 88 Ru, 90 Ru, 91 Ru, 92 Ru, 93 Ru and 94 Ru. 
The extended P + QQ model is confirmed to be useful 
in the heaviest N « Z nuclei. The results testify the en- 
hancement of the quadrupole correlations at the N = Z 
nucleus 88 Ru as compared with the other Ru isotopes. 

However, the disagreement between theory and exper- 
iment for 88 Ru cannot be disregarded. The slope of the 
J — to graph showing the moment of inertia and the 
durable increase of angular frequency are not sufficiently 
reproduced for 88 Ru with the QQ force strength com- 



monly fixed to all the Ru isotopes. The theoretical anal- 
ysis suggests a further enhancement of the quadrupole 
correlations, and recommends us to use a stronger QQ 
force for 88 Ru. We have tried to strengthen the p — n 
QQ interaction in the two ways so as to conserve and 
not to conserve the isospin of eigenstates, and also to 
strengthen all the p — p, n — n and p — n parts of the 
isoscalar QQ force. Within a small enhancement, how- 
ever, there is little to choose between them in the present 
calculations. Anyway, the present study indicates a spe- 
cial enhancement of the quadrupole correlations in the 
N = Z nucleus 88 Ru. This is consistent with the large 
deformation of 88 Ru in contrast to 90 Ru which is pre- 
dicted by the projected shell model calculation on the 
deformed basis 0, [l6| . 

The requirement of the enhanced QQ force for 88 Ru 
possibly means that the configuration space should be ex- 
tended in our spherical shell model. We have investigated 
the contribution of the 2d 5 / 2 orbit which is expected to 
mix with the lgg /2 orbit through the large Q matrix ele- 
ment. The truncated space (2pi/2, 139/2, 2d 5 / 2 ) can eas- 
ily reproduce the slope of J — u> graph observed in 88 Ru. 
The result suggests that the 2g?5/2 orbit contributes to 
the quadrupole correlations, which supports that 88 Ru is 
deformed. Contrary to this, the same calculation requires 
a much smaller contribution of the 2d 5 / 2 orbit to 90 Ru. 
It is, therefore, likely that 88 Ru is deformed while 90 Ru is 
not largely deformed as known from the observed J — uj 
graphs. This situation still demands different QQ force 
strengths for 88 Ru and 90 Ru in our spherical shell model 
calculation. There is not a selfconsistent way to deter- 
mine the QQ force strength. An additional constraint, 
for instance with respect to the Q moment value, is nec- 
essary for it. The condition is the same for the treatment 
on the deformed basis 0, 0] . The deformation should 
be selfconsistently determined there. 

Our model with a single set of parameters, however, 
is capable of describing the difference between 88 Ru and 
90 Ru. The calculations have presented a useful knowl- 
edge of the structure of Ru isotopes. The contrast fea- 
tures of 88 Ru and 90 Ru owe to the a-like (T = 0) 
2p — 2n correlations depending on the shell structure 
in the I39/2 subshell nuclei. In the 90 Ru isotope with 
one extra neutron pair which does not join in the a-like 
2p — 2n correlations, the extra neutron pair aligns easily 
to J = 9/2 + 7/2 = 8 (T = 1) in lg 0/2 - In contrast to 
this, the a-like (T = 0) 2p — 2n correlations hinder the 
single nucleon-pair alignment coupled to J = 8 (T = 1) 
till the simultaneous alignments of proton and neutron 
pairs at J — 2 x 8 (T = 0), in the N = Z even-even 
nucleus 88 Ru. 

The shell structure produces characteristic bands with 
opposite parities in 89 Ru. The following predictions are 
obtained for 89 Ru: The l/2~ state is the ground state. 
There are three characteristic bands. The negative par- 
ity band Nl on l/2~, which resembles the ground-state 
band of 88 Ru, shows backbending at 33/2~ caused by 
the simultaneous alignments of proton and neutron pairs 
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coupled to J — 16 in \g$/2- The J < 15/2 states of an- 
other negative parity band N2 on 3/2~ are the partners 
of the J < 17/2" states of the band Nl. The band N2 
shows a unique backbending at 19/2 - caused by the p — n 
pair alignment coupled to J — 9 (T = 0) in 1.99/2- The 
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